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Os '. Abstract 

We outline a path integral derivation of both the transverse force, the Berry 
' phase, and friction for a vortex from the microscopic fermionic superfluid 

theory. The derivation manifests transparently the mutual independence of 
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the Berry phase and dissipative terms in the effective vortex action. [] 



I. INTRODUCTION 

Vortices in superconductors and He3 superfluids are topological excitations. They de- 
termine the global properties such as the stability of the supercurrent carrying states, and 
have been under intensive theoretical and experimental studies since earlier sixties [|]]. The 
q . detailed microscopic Bardeen-Cooper-Schrieffer type (BCS) theory for fermionic superfluids 



is successful and well defined. The derivation of vortex dynamics is, however, non-trivial 
and less certain. Here we present a path integral derivation of vortex dynamics following the 
line of geometric methods |||3| . The advantage of the present method is that in the effective 
vortex Lagrangian the separation of geometric and dissipative contributions is natural. The 
important results here are that the transverse force agrees with the one obtained by the 
Berry phase method and by the total force- force correlation function method ||, and 
is insensitive to details. The friction is determined by the spectral function of the Hamil- 
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tonian, sensitive to details. The large transverse force has been found recently by a direct 
measurement ||, in consistent with the prediction in Ref. @,[|. 



II. DESCRIPTION OF THE BERRY PHASE AND DISSIPATION 

In this section we review the description of the Berry phase and dissipation for a particle 
within a model Hamiltonian, and show how those two quantities appear independently in 
the effective particle action. 

To study the quantum dissipative dynamics of a particle, one may specify the whole 
Hamiltonian, particle plus environment. In practice, what normally needed is the particle 
dynamics with the influence from the environment. The environment itself is left unobserved. 
Hence an accurate and efficient method to integrate out environmental degrees of freedoms is 
needed. This can be most readily done using the path integral method. In the present paper 
the approach and notation of Ref. || to this question will be followed. In that formulation, 
the dissipation is produced by an environment consisting of a set of harmonic oscillators. 
The corresponding total model Hamiltonian may have the following form: 

H = ^-[P + A(x)] 2 + [/(x)+E 

v j 

Here x = (x, y), P, {q^}, and {pj} are all two dimensional vectors. The vector potential 
A determined by V x A = B corresponds the transverse force — x x B, with B in the 
z-direction. The effect of the dissipative environment is specified by the spectral function 
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In the present paper, we shall assume the spectral function to have the following form 

J(uj) = rjuj s exp < 1 , (3) 

I u c ) 

with u c the cutoff frequency whenever needed. In accordance with Ref. |5|], s > 1 is the 
superohmic case, s = 1 the ohmic case, and < s < 1 the subohmic case. In the ohmic 
damping case, 77 is the friction coefficient in the usual Langevin equation. 
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In the imaginary time formulation the particle dynamics is described by the Euclidean 
action 



S 
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dr 
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where f3 = l/A^T is the inverse temperature. 

Since we are interested in the particle dynamics, keeping other degrees of freedoms of the 
environment is unnecessary. The integrations over {qj} will then be taken. The resulting 
effective action for the particle is MM 
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with the damping kernel F\\ as 
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It is evident that the transverse force corresponds the term ix • A in the (effective) particle 
action. It has not been influenced by the dissipative environment, due to the rotational and 
translational invariances of the coupling between the particle and the environment. This 
term gives the Berry phase, or the Aharonov-Bohm phase for a charged particle, if moving 
along a closed trajectory. The influence of the environment is contained in the non-local 
term. It should be pointed out that in the dissipative dynamics described by Eqs.(5,6) all 
one needs to know is the spectral function J{u). The details of the original environmental 
Hamiltonian have been suppressed. This implies that the dissipation can be produced by a 
fermionic environment, an observation important in our following derivation. 



III. EFFECTIVE VORTEX ACTION 



To find the effective vortex action, we begin with the standard BCS Lagrangian for s- 
wave pairing in the imaginary time representation. Here the unwanted fermionic degrees of 
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freedoms will be integrated out, instead of those bosonic ones in section 2. We will only 
consider neutral superffuids here, but the coupling to electromagnetic fields does no affect 
our main results. More detailed analysis will be published elsewhere. J7| The model BCS 
Lagrangian is 

-,2 
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]T i>t(x, t) (hd T -fi F - |- V 2 + V(x)^j if> a (x, r) 



(7) 



where describes electrons with spin o = (f, j), A*f the chemical potential determined by 
the electron number density, V(x) the impurity potential, and x = (x, y,z). A vortex at x v 
has been implicitly assumed. The partition function is 



Z = Jv{x v ,^,ip} x expj-i^ dr J d 3 xL BCS ^j , 
with (3 = 1/kgT , and d 3 x = dxdydz. Inserting the identity in the functional space, 

1 = Jv{A*,A}expl-^£ P drjd :i 
into Eq.(8) we have 



(8) 
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Z = i V{x v , if>\ij>, A*, A} x expl-^J ^ dr / tf.r x 
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Here the Hamiltonian is defined as 



1 r h P 
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dr / d 3 x\A\ 2 
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with H = —(h 2 /2m)\7 2 — fiF + V(x). Integrating out the electron fields and ip a first, then 
the auxiliary (pair) fields A under the meanfield approximation, one obtains the partition 
function for the vortex 



V{x v } exp 



S, 



(10) 



with the effective vortex action 



Seff - TrlnG- 1 - ±- [** dr [ d 3 x\A\ 2 , (11) 
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where Tr includes internal and space-time indices, and the Nambu-Gor'kov (NG) Green's 
function G defined by 

(hd T + H)G(x, t; x', t') = 6(t - t')5 3 (x - x'), (12) 

together with the BCS gap equation, or the self-consistent equation, 

A(x,r) = -g < V>|(x,t)V>t(z,t) > . (13) 

A special attention should be paid to the equal time limit of the NG Green's function. M 

We assume that the vortex is confined to move in a small regime around a point at xq, 
which allows a small parameter expansion in terms of the difference between the vortex 
position x v and xq. We look for the long time behavior of vortex dynamics under this small 
parameter expansion. For the meanfield value of the order parameter, this expansion is 
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A(x, r, x v ) = 1 1 + 5x v {t) ■ V + ^(Sx v (r) ■ V ) 2 J A (s, x ) . (14) 

Here 6x v = x v — Xq. In Eq.(14) we have used the fact that when x v = Xq the vortex is static. 
The effective vortex action to the same order is, after dropping a constant term, 

= __ Tr ( GoS ')2 + —J o dr J d 3 x x 5x v ■ V A* 6x v ■ V A , (15) 

with 



/ 5x v - V Ao ^ 
y Sx v ■ V A* j 



(16) 



Here Go is the NG Green's function with A(A*) replaced by A (Aq), and the gradient Vo 
is with respect to x . 

Before constructing the NG Green's function Go we consider the eigenfunctions of the 
Hamiltonian 7Y = W(A ,Aq). The stationary Schrodinger equation, the Bogoliubov-de 
Gennes (BdG) equation, is 



Ho^k(x) = E^ k (x) , (17) 



with 
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Since 7io is hermitian, all its eigenstates form a complete and orthonormal set. Eq.(17) have 
two useful properties: 



Hq^(x) = E V(x) => Ho^(x) = -E *(x) , 



(18) 



with 



and, since 5x v ■ Vq'Hq = for k ^ k', 
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J d 3 xVl(x)W k ,(x) = [E k , - E k )5x v ■ J d 3 xVl(x)V y k ,(x) 



(19) 



This equation implies that the two ways of specifying the vortex coordinate, through the 
trapping potential or through the order parameter, are equivalent. 
The NG Green's Go can be expressed as 



_i -io; n (r-r') 

Go(x,t;x',t') = J:—— -MxMW)- 
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Here uo n = nir/h(3, with n odd integers. Assuming the rotational symmetry after the impurity 
average, a straightforward calculation leads to the following effective vortex action 

1 rh/3 



S eff = - dr K \5x v (t)\ 2 + / dr'F||(r - r')\Sx v (r) - Sx v [ 
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with the spring constant in the effective potential, 

K = - [ d 3 x\VoA*Jx,x )\ 2 - (°° du^ 
g J Jo oj 
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the damping kernel, 
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and the transverse kernel, in the long time limit, in terms of the virtual transitions, 
F ± (t) = -d T - T ,5(r -r')E / ^ / rfV h (fk ~ fk>) 

k,k' 

or in terms of the contribution from each state, 

F±(r) = -d T - T ,6(r - r') W rf 3 x ftz • (/ fe Vo^(x) x V u k (x) 

k J 

-(1 - fk)V v* k (x) x V u fc (a;)) . 
To reach Eqs. (21-25), following two identities have also been used: 

P fk , 5 = o- 
-j3 (l-f k ), 5 = 0+ 
and 
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with the Fermi distribution function f k — 1/(1 + e^ Ek ), and the spectral function 



(27) 



J M = ?-E 5 (^- l^fe-^fe'l)l/fe-/fc'l x 
2 fc.fe' 
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In the following we discuss the implications of K , Fy , and Fj_ one by one, and show that 
Eq.(21) contains both the dissipative effect and the transverse force, or the Berry phase. 



IV. DISCUSSIONS 



For the purpose of getting the friction and the transverse force, the precise value of K 
in Eq.(22)is irrelevant. 
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The longitudinal kernel Fu contains all information on the vortex friction, revealed by 
the fact that Eqs. (21,23) are identical to the nonlocal action of Eqs.(5,6). The friction is 
determined by the low frequency behavior of the spectral function J(&), Eq.(28). 

If we classify the eigenstates of the BdG equation, Eq.(17), according to core(localized) 
and extended states, the dissipation comes from all three parts: core-to-core, core-to- 
extended (or extended-to-core), and extended-to-extended transition rates. Those rates are 
added up, because they are all positive, which differs from the situation for the transverse 
kernel. 

For a clean superffuid with no impurities, it is straightforward to find the damping 
corresponds to super-Ohmic cases. The dissipation is weak in this case. The main effect is 
to renormalize the vortex mass. j7| 

The presence of impurities mixes up all the clean limit eigenstates, and generates a quasi- 
continuous core energy spectrum after the impurity average. A perturbative calculation 
shows that it is indeed the case, and the spectral function becomes Ohmic. [^] This suggests 
that in the clean limit the dissipation is super-Ohmic, and turns into Ohmic in the dirty 
limit. 

In the long time limit there are two equivalent forms for the transverse kernel F±: The 
virtual transition expression, Eq.(24), and the individual state contribution, Eq.(25). Their 
equivalence has already been discussed in Ref. @ with the aid of conservation laws. Eq.(24) 
can be again casted into various equivalent forms because of the cancelations among those 
virtual transitions. For example, at zero temperature only two core states closest to the 
Fermi surface, one below '-' and one above '+', will be involved: 

F x = -d T „ T ,6(r - rO I d 3 x I d 3 x' &(/_ - /+) 

z ■ (tyl(x')V y+(x') x V V\.(x)V-(x)) • ( 29 ) 

In this case it is a statement of the spectral flow ||. One may also cast it into the forms of 
core-to-extended, or extended-to-extended transitions. 

In Eq.(25), the counting of individual state contributions is expressed as an area integral 
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of the momentum commutator. It can be expressed as the line integral of the momentum 
density far away from the core, 



with the momentum density j(x) = ih\im x ^ x > |(V — V')pi(x, x'), and the one-body density 
matrix 



Then F± corresponds the term, the Berry phase ih(p s /2) J dr 5x v -(Sx v xz)/2, in the effective 
vortex action Eq.(21). Here p s is the superfluid number area density. The immediate 
conclusion, as drawn in Ref. is that localized core states do not contribute, because both 
the total momentum density at the core Xq and the momentum density of core states far 
away from xq are zero. Thus we have two complete different ways to evaluate the transverse 
force: via virtual transitions at the core, or a counting of momentum density for each state 
far away from the core. This is similar to the evaluation of the quantum Hall conductance: 
the two methods, through either edge or bulk states, are equivalent. 

The insensitive to details such as impurities is most transparent from Eq.(25) or (31): 
In the one-body density matrix both the electron number density and the phase 6(x — x v ) 
are all insensitive to details, if the localization effect caused by impurities is negligible. Here 
the phase 9 is defined through the order parameter A(x,t,x v ) — > \A\e ig9( - x ~ Xv \ with q = ±1 
describing the vorticity along the z direction and 9{x) = arctan(?//x). 
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pi(x,x') =^2((f k ul(x)uk(x') - (1 - fk)vl(x)v k (x')) . 
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